NASA Technical Memorandum 85907 


N 8 4 - 2 7 7 8 b 


Thermal Response of Space Shuttle 
Wing During Reentry Heating 

Leslie Gong, William L. Ko, and Robert D. Quinn 


June 1984 


NASA 

National Aeronautics and 
Space Administration 





NASA Technical Memorandum 85907 


Thermal Response of Space Shuttle 
Wing During Reentry Heating 

Leslie Gong, William L. Ko, and Robert D. Quinn 

NASA Ames Research Center, Dryden Flight Research Facility, Edwards, California 93523 


1984 


iw\s/\ 

National Aeronautics and 
Space Administration 
Ames Research Center 
Dryden Flight Research Facility 
Edwards, California 93523 





THERMAL RESPONSE OF SPACE SHUTTLE WING DURING REENTRY HEATING 


Leslie Gong,* William L» Ko, * * and Robert D. Quinnt 
NASA Ames Research Center 
Dryden Flight Research Facility 
Edwards, California 93523 


Abstract 

A structural performance and resizing (SPAR) 
finite element thermal analysis computer program 
was used in the heat transfer analysis of the 
space shuttle or biter that was subjected to 
reentry aerodynamic heatings . One wing segment 
of the right wing (WS 240) and the whole left wing 
were selected for the thermal analysis. Results 
showed that the predicted thermal protection sys^ 
tern (TPS) temperatures were in good agreement with 
the space transportation system, trajectory 5 
( STS -5 ) flight-measured temperatures. In addition, 
calculated aluminum structural temperatures were 
in fairly good agreement with flight data up to 
the point of touchdown. Results also showed that 
the internal free convection had a considerable 
effect on the change of structural temperatures 
after touchdown. 

Introduction 

The space shuttle orbiter reenters the earth’s 
atmosphere at an altitude of about 121,920 m 
(400,000 ft) and at extremely high velocity 
(nearly Mach 25 at the time of reentry). In order 
to protect the shuttle structure from severe 
reentry aerodynamic heating, the entire shuttle 
structure is covered with a thermal protection 
system (TPS). The regions of the shuttle surfaces 
receiving lower heating rates (for example, wing 
upper surfaces, fuselage side walls, and bay doors) 
are covered with highly flexible felt reusable 
surface insulation (FRSI), and the regions exposed 
to higher heating rates (for example, wing and 
fuselage lower surfaces) are covered with TPS 
tiles. A layer of highly flexible strain isola- 
tion pad (SIP) is sandwiched between the TPS tiles 
and the aluminum skin for absorption of the strain 
incompatibility between the brittle TPS tiles and 
the aluminum skin. Overheating of the aluminum 
structure may cause thermal creep resulting in the 
loss of structural integrity required for subse- 
quent flights. Although the SIP layer may absorb 
the aluminum skin thermal buckling effect on the 
TPS tiles to some extent, the excess thermal buck- 
ling of the aluminum skin because of overheating 
could cause a debonding of TPS tiles resulting in 
partial or total loss of the TPS protecting func- 
tions. In former space transportation system (STS) 
flights the performance of the TPS was so good 
that the shuttle structural temperatures during 
reentry were kept well below the design limit tem- 
perature of 176° C (350° F) , and thus the above 
mentioned concerns were practically eliminated. 
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Since the shuttle is designed to be used as 
many as one hundred times, it is essential to 
understand the mechanical performance (that is, 
structural stress levels) under the reentry aero- 
dynamic and thermal loadings so that confidence in 
the shuttle structural integrity can be established. 
The flight load data obtained from onboard strain 
gage measurements contain both the mechanical and 
thermal components. Unfortunately these two com- 
ponents cannot be separated experimentally. In 
order to obtain mechanical stress, thermal stresses 
must be removed from the strain gage measured 
stresses. This can be done analytically by first 
calculating the thermal stresses and then removing 
the thermal stresses from the strain gage measured 
stresses to give the true mechanical stresses. To 
calculate the thermal stresses, the structural tem- 
perature distribution must be known. Since the 
number of onboard thermocouples are insufficient 
to accurately record the structural temperature 
distribution, heat transfer analysis must be per- 
formed to calculate accurate structural tempera- 
ture distribution which will be used as input to a 
structural model for thermal stress calculations. 

The preliminary heat transfer analyses of typical 
wing cross sections were reported in Refs. 1 , 2, 
and 3. This paper presents results of an exten- 
sion of those analyses and discusses finite ele- 
ment heat transfer analysis of a wing segment of 
the right shuttle wing and the whole left shuttle 
wing. The results are compared with the STS-5 
temperature data (the most complete set of STS 
flight data). 

Nomenclature 

A^ surface area of radiation exchange ele- 

ment i 

C capacitance matrix 

Fij radiation view factor, or fraction of 

radiant heat leaving radiation 
element i incident on radiation 
element j 

FRSI felt reusable surface insulation 

HRS I high temperature reusable surface 

insulation 

JLOC joint location (grid point) 

K k conduction matrix 

K r radiation matrix 

K21 two -node conduction element 

K31 three -node conduction element 

K41 four-node conduction element 



K61 

K81 

LRSI 


Q 

R 

RTV 


six-node three-dimensional conduction 
element 

eight-node three-dimensional conduction 
element 

low temperature reusable surface 
insulation 

source heating load vector 
radiation load vector 
room temperature vulcanized 


of 121,920 m {400,000 ft) „ The nominal (or 
designed) trajectory is shown with solid curves, 
and data points were obtained from the space trans- 
portation system, trajectory 5 (STS- 5) flight. 

The trajectories for STS-1 through STS- 4 are simi- 
lar to the STS -5 flight trajectory and, therefore, 
are not shown* 

The calculation of the reentry aerodynamic 
heating was based on the actual STS-5 flight tra- 
jectory. The STS-5 flight was chosen because it 
provided the most complete set of flight data com- 
pared with other STS flights in which some of the 
flight data were lost. 


R31 


three-node radiation element 


Description of Wing Structures 
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SIP 

SPAR 

STS 

T 

THEOSKIN 
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t 

WS 
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four-node radiation element 

strain isolation pad 

structural performance and resizing 

space transportation system 

absolute temperature 

theoretical skin heating 

thermal protection system 

time 

wing segment 
longitudinal station 
wing station 
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Description of Problem 


Figure 1 shows the location of the wing seg- 
ment WS 240 of the right shuttle wing selected 
for the heat transfer analysis. WS 240 is bounded 
by yg 240 and yg 254 (Fig. 1). The leading edge 
portion of the wing and the elevon were not 
included in the analysis. Figure 2 shows the left 
shuttle wing that was analyzed. The region of the 
whole wing was bounded by longitudinal stations 
xo 1040 and xq 1365 (that is, elevon, glove, and 
leading edge regions were not included in the 
analysis) . * 

The reentry trajectory for the shuttle is 
shown in Fig. 3. The time was counted from the 
start of the reentry that occurred at an altitude 


*The existing structural model for the thermal 
stress calculations does not include the secondary 
load carrying structures, elevon, leading edge 
regions of the wing, and the glove. Therefore, in 
the thermal analysis, only the major load carrying 
structures were modeled. 


As shown in Fig. 2, the ribs of the wing are 
fabricated with vertical and diagonal trusses. 

All of the spar webs and the wheel well vertical 
walls are made of corrugated aluminum skins . Both 
the lower and upper wing skins lying between 
Xq 1191 and x Q 1365 are made of hat-stringer- 

reinforced aluminum panels. The lower and upper 
wing skins lying between xg 1040 and x 0 1191, and 
the leading edge beam assembly are made of alumi- 
num honeycomb -core -reinforced sandwich panels. 

The landing gear door is made of double -walled, 
hat-stringer-reinforced aluminum panels separated 
by aluminum stringers. The entire lower wing is 
covered with high temperature reusable surface 
insulation (HRS I) tiles, with a layer of SIP lying 
between the wing skin and the HRSI for absorption 
of thermal strain incompatibility between the wing 
skin and the HRSI. Most of the upper skin near 
the leading edge region is covered with low tem- 
perature reusable surface insulation (LRSI) tiles. 

A SIP layer lies under the LRSI. The rest of the 
upper wing skin, which is subject to low heating, 
is covered with FRSI under which there is no SIP 
layer. The room temperature vulcanized (RTV) 
rubber bonding agent was used in bonding the TPS 
and SIP to the substructures . Some of the gaps 
between the TPS tiles in the high temperature areas 
were filled with ceramic coated alumina mat (gap 
fillers) to prevent hot gases from coming in con- 
tact with the substructure at the bottom of each 
gap. Figure 4 shows the structural details of the 
wing segment WS 240 of the right shuttle wing. 

Thermal Modelings 
Structural Simplifications 


Because of the complex nature of the shuttle 
structure, some structural simplifications were 
required before setting up thermal models so that 
the analysis would be manageable using existing 
computers. Excessively detailed models could lead 
to tedious radiation-view-factor computations, in 
which the gain in solution accuracies may not be 
as much compared with the solution obtained from 
simpler, yet reasonably detailed models. In order 
to examine the adequacy of representing the hat- 
stringer -reinfo reed skins with smooth skins of 
uniform equivalent thicknesses, the conventional 
finite difference method was used in the two- 
dimensional heat transfer analysis of a single 
hat-stringer. As shown in Fig. 5, the hat-stringer 
that was analyzed was located at the lower skin of 
WS 240, bay 3. The upper skin, the spar webs, and 
the lower skin — excluding the hat stringer — 
were assumed to have uniform thicknesses. In the 
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analysis, all types of radiation heat exchanges 
were considered: (1) external radiation from the 
TPS surface into space; (2) radiation exchanges 
between the hat-stringer outer surface and the 
inner surface of the bay; and (3) internal radia- 
tion inside the hat stringer « The heat input was 
based on that used on Mission 3. The calculated 
temperatures in Fig, 6 show the peak temperature 
difference between points A and B to be about 
14. 44° C (26° F) , The temperature curves for the 
case where there was no internal radiation within 
the hat stringer is almost identical to correspond- 
ing temperature curves in which the hat -stringer 
internal radiation was considered. This suggests 
that the effect of the hat-stringer internal radia- 
tion is negligible. Because the temperature dif- 
ferences between points A and B of the hat stringer 
were small, in the finite element modeling, the 
hat stringer, corrugated spar webs, wheel well 
walls, and honeycomb-core sandwich panels were 
represented by smooth solid panels of equivalent 
thicknesses having isotropic thermal conductivities 

Finite Element Modelings 

A structural performance and resizing (SPAR) 
finite element thermal analysis computer program 
(Ref. 4) was used for the reentry heat transfer 
analysis of the space shuttle orbiter. 

Wing Segments . The SPAR finite element model 
setup for the wing segment WS 240 is shown in 
Fig. 7. The model is three-dimensional and has 
920 joint locations (JLOCs), The wing skins, 
spar webs, rib cap shear webs, RTV layers (lying 
on both sides of the SIP), and the TPS surface 
coatings were modeled with K41 elements (four -node 
heat conduction elements). The spar caps, rib 
caps, and the rib trusses were modeled with K21 
elements (two-node heat conduction elements). The 
TPS was modeled in 10 layers on the lower surface, 
and 3 to 4 layers on the upper surface using K81 
elements (eight-node three-dimensional heat con- 
duction elements) and K61 elements (six-node three- 
dimensional heat conduction elements). The K61 
elements were used only in the region where the 
modeled TPS layers changed from four to three lay- 
ers on the bay 1 upper surface (see Fig. 7). The 
SIP was modeled with only one layer of K81 ele- 
ments . The aerodynamic surfaces were modeled 
using one layer of K41 elements of unit thickness 
to provide source heat generation. The external 
and the internal radiations were handled by attach- 
ing a layer of R41 elements (four-node radiation 
elements) to the radiation surfaces. The radia- 
tion into space was modeled by one R41 element. 

No radiation elements were attached to the sur- 
faces of the rib cap shear webs and the rib trusses 
because of small exposed areas. The front and rear 
ends of the model were totally insulated. 

Because of the existence of the gaps between 
the TPS tiles HRS I and IRS I, the heat flow through 
the TPS tiles was restricted only in the tile 
thickness direction for the thermal models men- 
tioned above. In the analysis, two thicknesses 
were used for the HRSI and LRSI tiles: 100-percent 
and 80-percent TPS thicknesses. The purpose of 
using the effective thickness of 80 percent of the 
original TPS thickness was to account for the gap 
heating between the TPS tiles. The effect of 


internal natural convective heat transfer was not 
considered.** The effect of the forced convective 
cooling on the structural temperatures prior to 
touchdown was found to be negligible and, there- 
fore, was not considered for the wing thermal 
models except the WS 240, 

Whole Wing . The SPAR thermal model (called 
WING model) set up for the whole left shuttle wing 
is shown in Fig. 8, This wing model was patterned 
after existing structural models. The entire WING 
SPAR thermal model has 2289 joint locations. 
Because of the exsistence of gaps between the TPS 
tiles, HRSI and LRSI, and because of the thinness 
of both FRSI and SIP layers, the heat flow through 
TPS and SIP was assumed to be one -dimensional in 
thickness direction. Thus the TPS for the lower 
surface was modeled in five layers with K21 ele- 
ments, and the TPS for the upper surface was mod- 
eled in three layers with K21 elements. Using K21 
elements was found to be more efficient (with 
respect to computation time) than using K81 ele- 
ments and restricting the heat flow to one direc- 
tion. The SIP layer was modeled with only one 
layer of K21 elements. The rib trusses and the 
rib and spar caps were modeled with K21 elements. 
All of the aluminum skins, the corrugated spar 
webs, the wing root wall (fuselage wall), and the 
wheel well walls were modeled with K31 and K41 
elements. The aerodynamic surfaces were modeled 
with one layer of K41 elements of unit thickness 
for source heat generation. The internal and 
external radiation effects were handled by attach- 
ing R31 (three-node radiation elements) and R41 
elements at the radiation surfaces of the aluminum 
skins and the TPS. The wing root vertical wall 
and the vertical boundary walls of the wing were 
completely insulated. The radiation into space 
was modeled with only one R41 element which was 
kept at a constant temperature of 26.7° C (80° F) , 
The effect of natural convection inside the wing 
was not considered. 

There were three purposes for setting up the 
WING SPAR thermal model for the whole wing. The 
first was to generate the temperature distribution 
required for input to the existing whole wing SPAR 
structural model for the calculations of thermal 
stresses. The second was to compare the predic- 
tions with those based on wing segment models. 

The third purpose was to examine the adequacy of 
the wing model because it was recognized that the 
WING model was not fine enough to define tempera- 
ture gradients near spars and ribs . The reason 
for using coarser elements for the whole wing 
model compared with those used for the wing seg- 
ment model was to match the existing structural 
model for the whole wing and to limit the effort 
in the computations of view factors for the three- 
dimensional radiation exchanges inside the wing. 

Aerodynamic Heating 


The external heat inputs to the thermal models 
were computed by a NASA computer program called 
theoretical skin heating (THEOSKIN) which uses the 
velocity-altitude-angle-of -attack time history of 


**At the time of this writing, the capability of 
handling two-dimensional free convection is being 
introduced to the SPAR program. 
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the flight-measured STS-5 trajectory given in 
Fig, 3* The THEOSKIN program solves the one- 
dimensional thin skin heating equation and com- 
putes time histories of surface temperatures, 
heating rates, heat transfer coefficients, and 
skin friction. The thermodynamic and transport 
properties of air used in this analysis are given 
in Ref® 5. 

Representative heating rates for WS 240 are 
presented in Fig, 9. The heating rates for the 
lower surface were computed assuming laminar flow 
up to 1160 sec and turbulent flow from 1160 sec 
until the end of the flight. The laminar heat 
transfer was computed by relating heat transfer to 
a skin-f riction equation through a modified Rey- 
nolds analogy. In this analysis, the Blasius 
incompressible skin friction equation was related 
to heat transfer by the Prandtl number to the neg- 
ative six-tenth power.® Compressibility effects 
were accounted for by using Eckerts* reference 
enthalpy transformation. 7 Details of this method 
for calculating heat transfer at hypersonic speeds 
are given in Ref. 8. The turbulent heat transfer 
was computed by a similar procedure except that the 
van Driest transformation was used to account for 
compressibility, and the Reynolds analogy factor 
was assumed to be a constant value of l.l.^' 1 ® 

The boundary layer flow on the upper surface 
of bay 1 at WS 240 was assumed to be attached. 

The rest of the upper surface of the wing was 
assumed to be in a region of separated flow. The 
heat transfer for the attached flow areas was com- 
puted using the same heat transfer codes that were 
used to calculate the lower surface heating. To 
calculate the heating rates for the separated flow 
areas on the upper surface, the heat transfer codes 
were empirically modified. The empirical correc- 
tions that were used were determined from compari- 
sons with previous measured flight data. 

The heating rates for the whole wing were gen- 
erated from calculations made for WS 240 and other 
wing segments through linear interpolations and 
extrapolations . 


Radiation Exchange 


For both the external and the internal thermal 
radiation exchanges, all of the view factors were 
calculated from the equation 


A i F ij “ A j F ji 


( 1 ) 


where A^ is the surface area of radiation exchange 
element i, and Fjj_ is the view factor, defined as 
the fraction of radiant heat leaving element j 
incident on element i. 11 In the calculation of 
view factors for the external radiation exchanges 
(considering that element i represents the space 
element and element j any radiation exchange ele- 
ment on the wing or fuselage surface), Fj^ was 
taken to be unity; therefore, Fij = Aj/Aj^ accord- 
ing to Eq. (1). Values of emissivity arid 
reflectivity used to compute radiant heat fluxes 
were: 


Surface 


Emissivity Reflectivity 


Windward 0.850 
Leeward 0.800 
Internal structure 0.667 
Space 1.000 


0.150 

0.200 

0.333 

0 


The initial temperature distribution used in the 
analysis was obtained from the actual flight data. 
In thermal modeling, the majority of the time was 
consumed in the computations of view factors. 

Transient Thermal Solutions 

The SPAR thermal analysis finite element com- 
puter program was used in the calculation of tem- 
perature time histories at all joint locations of 
the thermal models. The SPAR program used the 
approach described below to obtain transient ther- 
mal solutions. 

The transient heat transfer matrix equation 

(% + K r )T + CT=Q+R (2) 


where 

= conduction matrix 
K r = radiation matrix 

T = temperature 
C = capacitance matrix 


Q =s source load vector 
R = radiation load vector 
t = time 


was integrated by assuming that the temperature 
vetor T i + 1 at time step t^ +1 can be expressed as 

T i + 1 = Ti + TjAt + Jj-Td-At 2 + At3 + ••• O) 

where 


[”] 


a 2 [ 1 

at 2 


f '] 


a 3 [ ] 

at 3 


and T^ is the temperature vector at time step t^ 
and At is the time increment. The vector T^ is 
determined directly from Eq. (2); that is. 


Ti = -C -1 (K k + K r )Ti + C _1 (e + R) (4) 

Higher -order derivatives are obtained by differen- 
tiating Eq. (2) according to the assumption that 

(1) material properties are constant over At, 

(2) Q varies linearly with time, and ( 3) R is con- 
stant over At. For example. 


Ti = -C -1 (K k + 4K r )Ti + C _1 Q (5) 

Ti = -C -1 (K k + 4K r )T i + 4K r Ti (6) 

In the present computations, the Taylor series 
expansion (Eq. (3)) was terminated after the third 
term. The pressure dependency of the TPS and SIP 
thermal properties was converted into time depen- 
dency based on the trajectory given in Fig. 3. 

Time -dependent properties were averaged over "time 
intervals" (RESET TIME) which were taken to be 2 
or 25 sec. Temperature -dependent properties were 
evaluated at the temperatures computed at the 
beginning of each time interval. Q, Q, and R were 
computed every 2 sec. 
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Results 

TPS Surface Temperatures 

Figures 10 to 13 show the locations and num- 
bers of thermocouples at the TPS surfaces and the 
aluminum skins » 

In Fig* 14 the STS-5 flight-measured TPS sur- 
face temperatures are compared with the calculated 
values at the wing cross section Yq 240. The 
agreements are quite good, indicating that the 
calculations of the aerodynamic heatings were 
accurate. In a narrow time range of t ~ 1500 to 
2000 sec (immediately before and after the touch- 
down time), the flight data, especially for the 
lower TPS surfaces, showed temperatures lower than 
the calculated temperatures. This discrepancy 
could be caused by insufficient forced convective 
coolings in the heat input calculations and the 
neglect of internal free convection resulting pri- 
marily from the entering of cool air into the 
shuttle. t Results showed that use of the forced 
convective cooling near the touchdown time (nega- 
tive heating) resulted in excess computation time 
because of the change of sign in heat input, and 
that the effect of using negative heating on the 
structural temperatures was almost pictorically 
inconspicuous. Therefore, negative heating was 
not included in the thermal analyses of the whole 
wing. The flight data for the lower TPS surface 
of WS 240, bay 3 contain void data as indicated, 
because the lower limits of the thermocouple tem- 
perature readouts were set too high. 

Structural Temperatures 

Figure 1 5 shows the calculated aluminum skin 
temperatures at typical points of WS 240, compared 
with the STS-5 flight data. For the wing segment, 
the temperature calculations made were based on 
80-percent TPS thickness (solid lines) and 100- 
percent TPS thickness (dashed lines). For the 
whole wing, 100-percent TPS thickness was used in 
the temperature calculations. For WS 240 (Fig. 15), 
the measured lower skin temperature data follows 
fairly close to the calculated temperature curves 
based on 100-percent TPS thickness up to almost 
touchdown time; after that the flight data contin- 
ued to show lower temperatures. The structural 
temperatures calculated from the wing segment 
model WS 240 (100-percent TPS thickness) and from 
the WING model are quite close. The marked dis- 
crepancies between the calculated and the measured 
lower skin temperature after touchdown could be 
caused by the effects of free convective cooling 
and the entering of cool air into the wing. The 
effect of free convections was not considered in 
the analysis because at the time of the analysis, 
the SPAR program was still incapable of handling 
free convective heat transfer. For the upper skin 
of yo 240, the measured and the calculated tempera- 
tures (based on 100-percent TPS thickness) compare 
reasonably well even after the touchdown. The 
agreement was especially good for the bay 1 upper 
skin. 


tAir enters the interior of the shuttle or biter 
at an altitude of 30,048 m (100,000 ft) (about 
1400 sec from reentry). 


Figure 16 shows the chordwise distributions of 
aluminum skin temperatures for the reentry time 
of 1600 sec for WS 240. The "scalloped" shape 
reflects the drop in temperature at the heat sinks 
(that is, spar caps). The scalloped distributions 
of the skin temperatures are a direct correlation 
to the degree of thermal stress buildups in the 
structure. The STS-5 flight data correlate rea- 
sonably well with the predictions except for a few 
data points. 

Conclusions 

The structural performance and resizing (SPAR) 
finite element computer program was used in the 
reentry heat transfer analysis of one wing seg- 
ment and the whole left wing of the space shuttle 
orbiter. The SPAR finite element thermal models 
were set up in three dimensions. The results of 
the thermal analysis showed that; 

1 . The predicted thermal protection system 
(TPS) surface temperatures agree quite well with 
the flight measured data. 

2. The measured and predicted structural 
temperatures correlate fairly well prior to touch- 
down. This implies that the SPAR thermal models 
set up for both the wing and the wing segment were 
adequate . 

3* The internal free convection has con- 
siderable effect on the structural temperatures 
after touchdown and cannot be ignored. 

4. The effect of the internal radiation was 
found to be significant and cannot be ignored even 
at the relatively low structural temperatures. 

The view factor computations for the internal 
radiation resulted in a major task in thermal mod-^ 
eling; therefore, capability of view factor com- 
putations internally in the SPAR program is highly 
recommended • 

5. Although the finite elements used in the 
SPAR WING model are relatively coarse, the thermal 
model gives fairly good structural temperature pre- 
dictions. Thus the structural temperature distri^ 
bution obtained from the WING model could be used 
with confidence for input to the whole wing SPAR 
structural model for thermal stress calculations. 

6. The SPAR WING thermal model gives only 
coarse structural temperature distribution for the 
overall estimation of thermal stresses in the 
whole wing. For the detailed local thermal stress 
calculations, the wing segment models give better 
solutions • 
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Fig. 6 Temperature time histories at points in hat-stringer mission 3 
heating • 
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Fig. 7 Three-dimensional SPAR finite element thermal model for WS 240 
(920 joint locations ) • 


Fig. 5 Space shuttle wing SPAR thermal model (2289 joint locations ) 
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Fig. 10 Wing lower TPS surface joint locations and thermocouple 
locations . 
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Fig. 11 Wing upper TPS surface joint locations and thermocouple 
locations . 






































Fig , 12 Wing lower skin joint locations and thermocouple locations . 



Fig* 13 Wing upper skin joint locations and thermocouple locations • 
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Fig. 25 JVS 240 structural temperature time histories, STS- 5 flight 
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